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COMPUTATION OF INFRARED COOLING RATES 
IN THE WATER VAPOR BANDS 

ABSTRACT 

A fast hut accurate method for calculating the infrared radiative terms due 
to water vapor has been developed. It makes use of the far wing approximation 
♦o scale transmission along an inhomogeneous path to an equivalent homogene- 
ous path. Rather than using standard conditions for scaling, the reference tem- 
peratures and pressures are chosen in this study to correspond to the regions 
win re cooling is most significant. This greatly increased the accuracy of the new 
method. Compared to line by line calculations, the new method has errors up to 
4 7c of the maximum cooling rate, while a commonly used method based upon 
the Goody band model (Rodgers and Walshaw, 1966) introduces errors up to 
1 1%. The effect of temperature dependence of transmittance has also been 
evaluated; the cooling rate errors range up to 1 1% when the temperature de- 
pendence is ignored. In addition to being more accurate, the new method is 
much faster than those based upon the Goody band model. 
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COMPUTATION ()l INI KAMI I) COOLING RAILS 
IN Till WATER VAPOR BANDS 

I. Introduction 

Absorption, emission, and transmission of thermal infrared radiation by water vapor 
have a very strong effect on the exchange of radiation energy between atmospheric layers 
and on the loss of radiation energy to space. Because of the influence of these radiative 
processes on the thermodynamic structure of the atmosphere, it is important both for 
modelling the general circulation and for climate modelling to accurately compute the terms 
involved. 

The most accurate methods for computing the radiative terms in a molecular atmosphere 
involve a detailed “line by line" calculation of the absorption coefficient versus wave number. 
In recent years, a fairly complete set of line parameters for the important atmospheric gases 
has been compiled and made available by McClatchcy et al. ( 1973). Although the underlying 
physics is well understood, the accurate line by line calculation of the radiative terms is a 
tedious, time consuming process which, for an atmospheric model, must be repeated at each 
grid point, every time step. Because of computing time limitations these methods are quite 
impractical for general circulation and climate models. Therefore, the band model and the 
emissivity method which can avoid the line by line detail are commonly used. Lor example. 
(iISS GC'M (Somerville et al.. 1974) partly makes use of the Goody band model to calculate 
transmit lances, and NC'AR GCM uses Hssascr and Culbertson’s ( I960} absorptivity approach 
with transmittances computed at a fixed temperature (Sasamori, 1968). The diffuse approxi- 
mation is used to compute diffuse transniittances for both approaches. 


In (his study wc develop a new method for calculating the radiative terms involving 
water vapor specifically, the atmospheric cooling rates which is both fast and accurate. 
The approach is based upon ( 1 ) the actual distribution of line strength and spacing as de- 
rived from the line parameter data; and (2) a scaling approximation based upon the behavior 
in the far wings of absorption lines, which dominate the spectrum. For application to cool- 
ing rate calculations in atmospheric models, the method can be made extremely fast by using 
precomputed tables. 

Tlte far wing approximation is developed in Sec. 2. Combined with the k-distribution 
method, it is applied in Sec. 3 to derive expressions for fluxes and cooling rates. In Sec. 4 
we compare cooling rate calculations based upon the new method, the Rodgers and Walshaw 
method, and the line by line method the last assumed to be the standard; we also test the 
validity of other approximations. Concluding remarks appear in Sec. 5. 



2. The far-wing approximation 

The computation of radiative fluxes and the atmospheric cooling rates involves the 
calculation of absorption due to thousands of molecular lines. The parameters that determine 
the absorption coefficient at any wave number v are the strength. S. half width, a. position 
v ' , and shape for each line in the spectrum. The dependence of line strength on temperature 
varies with each line, while the half width usually has a fixed relationship with temperature 
and pressure. In computing fluxes over an inhomogeneous path, both line strength and line 
width change along the path, and these effects must be taken into consideration. 

In the troposphere and lower stratosphere, the broadening of lines is primarily due 10 
molecular collisions and therefore the line shape is adequately represented by the I orentz 
function. In that case, the absorption coefficient at wave number v, pressure p. and tempera- 
ture T is given by 


k„(p.T) 


I 


S,(T) o ( (p.T) 

S n — • 

n (y- v oi y +a 1 “(p,T) 


(II 


where the summation is over all absorption lines, denoted by the subscript i. the hall width is 

/ p\/T f Vi 

given by a(p,1) = a <p ( .1 r J .and the subscript r denotes reference values. 

For any one line, the absorption coefficient increases with decreasing pressure for 
|r> - n 0 | < ur and decreases for | v - 1 > 0 | > a. In particular, it is proportional to pressure, p. 
in tin far wing of a line and to I Ip at the line center, l or atmospheric conditions, the half 
widths of the absorption lines are small compared to the mean line spacing and. therefore, 
the w,ngs dominate the spectrum. Furthermore, at those atmospheric levels where water 
vapor cooling is significant, absorption quickly saturates near the centers of most of the 
water vapor lines; therefore, errors in the absorption coefficient in the vicinity of line centers 
will have a small effect on transmittance. These arguments suggest that one might 
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approximate the pressure and temperature dependence of the absorption coefficient by its 
behavior in the wings of the lines. At one extreme, where | v - | » «j, the absorption 

coefficient can be written as 


where 


k„(p,T) « ky(Pr*T,) ” R„(T.T r ) , 


R„(T.T r ) 
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Figs, la and b show curves of R^(T, T r ) for specified values of T r and T * T f ± 40 K It 
can be seen that the absorption coefficient is less sensitive to temperature in the band center 
region ( 10-340 cm* 1 ) than in the band wing region (340-580 cm - ’). Although the variation 
of R„(T, T f ) with wave number is large, it is still small compared to the variation of k^. If we 
replace R,,(T. 1 f ) by its mean value, R(T, T f ). over a wide spectral interval, Eq. (2) becomes 

k„(p,T) ■ k,(p r ,T r )^R(T,T r ). (4) 

The use of (4) to represent the temperature and pressure dependence of k^tp, T) will 
be called the “far wing approximation.” The advantage of this approximation is the separa- 
tion of v from the variables p and T. As will be seen later, it leads to a great simplification 
in the computation of the radiative terms. 

The errors introduced by the far-wing approximation will not be large because of the 
rapid decrease of water vapor mixing ratio with height. This causes the vertical profile of 
cooling rate due to any spectral region to have a sharp peak, thus spanning relatively narrow 
ranges of temperature and pressure. Within those narrow ranges. Eq. (4) would be expected 
to be fairly accurate. 
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The reference parameter* p r and I' r could be fixed for the entire water vapor spactrum 
or one can divide the spectrum into several regions, witli p f and I, chosen in each region to 
correspond to conditions at the level of maximum influence on cooling rates. We found 
that one can obtain very pood accuracy by dividing the water vapor spectrum (including 
boll* :he 6.3 /am and the rotational bands) into only two iy/n s of spectral regions ( 1 1 the 
vicinity of the hand centers, where absorption is strong, and i 2) the hand winn\. where ab- 
sorption is moderate. The assignment of spectral regions to band w.ngs and band centers is 
shown in Table I. The reference parameters we»e chosen to be p f = 275 mb. I f = 225 k 
for the band centers and p r = 550 mb. T f = 256 K for the band wings, they correspond to 
points on the temperature versus pressure curve for a mid-latitude standard atmosphere 
(see Fig. 2a). where each spectral legion type contributes most significantly to the cooling 
rate. Also tabulated in Table I are the sallies R at T f ± 40 K for the two band regions. As 
is well known, the table shows that the band wings are much more sensitive to temperature 
than the band centers. Rather than calculate R(T. T r ) at a large number of temperatures, 
we took advantage of its expected smooth behavior and used a polynomial fit to the values 
at T = T f - 40. T = T f ( lor which R = 1 1 and T = I , + 40. 

The computational advantage of the far wing approximation will become obvious in the 


next section when calculating the transmittance through an atmospheric layer averaged over 
a wide spectral interval. Meanwhile, the monochromatic transmittance between pressure 
levels pj and pj. for radiation in .• direction forming the angle cos" 1 n with respect to the 


zenith axis, is given by 


(VP,. P2>) 


-- u„(p,.p 2 ) 


where u,, <p ( . p-, i is the optical thickness of the layer, defined by 
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(6) 


k <p. T ( p>> q(p> dp. 


<1 is the specific humidity and g is the gravitational acceleration. 

In the far wing approximation, the optical thickness becomes 


u„IP|.Ps) * k„(p r .T r ) w(p,.p 2 ) 


where 


r Pj | 

w(Pj. p 2 ) = / -- — K ( Tip). T, )q (pi dp (8) 

Jp t * Pr 

is the amount of water vapor in the layer, scaled by the factor (p/p, ) R [T(p). T, J . The 
scaling factor is a function not only oi the ambient temperature and pressure but also ot the 
selected reference conditions, p, and T ( . This constitutes a one-parameter scaling approxi- 
mation. It differs in a rather significant way from the more commonly used scaled water 
vapor amount (Hlsasser and Culbertson. I960; Manabe and Moller, 1961. I acis and Hansen. 
1974; Kamanathan. 1976; Wang. 1976) 


— f’^ 


r, • 


dtp) dp. 


where 0 < m < I and n > 0. F:q. (9) implies that absorption decreases with increasing tem- 
perature. quite the contrary of what actually happens, as revealed in F ig. 1 and I able 1 Of 
the two factors in (3) which determine the temperature dependence of the absorption 
coefficient, the second one is effectively ignored by the use of (9). and that turns out to be 
the dominant factor. Another important difference between w( p ( . p, ) and w' i p ( . p, ) is 
the choice of reference levels: in utilizing v '(p. . p, ), the reference conditions were always 
chosen at I atmosphere, thus introducing a large error in *he middle and upper troposphere, 
where cooling due to water vapor is large < as will be seen below ) By judicious choice of 
reference values, the error is significantly reduced at the important levels. 



To test the accuracy of the far wing approximation, transmit tances were calculated 
monochromatically at .01 cm' 1 intervals for/i ■ I between typical atmospheric levels for 
the mid-latitude atmosphere shown in Figs. 2a and 2b. Fig. 3 shows the transmittances be- 
tween 100 and 250 mb (Fig. 3a) and between 250 and 400 mb (Fig. 3b) plotted against 
wave number in the range 320 to 340 cm' 1 (near band center), and Fig. 4 shows transmit- 
tance between 300 and <>00 mb for wave numbers 4(>0 to 480 cm" 1 ( band wing region). In 
each gTaph. the solid curve is an “exact” calculation based upon the Voigt profile, taking 
into account its pressure and temperature variation along the path. The dashed curve 
represents the “far wing approximation." in which the absorption coefficient is computed 
according to (4). As might be expected, the far wing approximation is better for those 
layers close to the reference conditions. 

Mean tran. mittance over 20 cm" 1 wave number intervals for various atmospheric layers 
is shown in Table 2. It is seen that even far from the reference conditions the error is ^ .01. 
Over wider spectral intervals the errors are smaller. 
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3. I 'luxes and cooling rales 


In a non-scattering atmosphere with temperature profile I (pi, the monochromatic 
upward and downward fluxes of thermal radiation due to water vapor at pressure level p 
arc given by 


J r p % 

B„(T(p')) 

D 


d f I Uj,( p, p'|) 

— dp- 

dp 


(Ida) 


J f 1 ’ , v df(u,,(p, p'|| 

B,(T<p'>) dp' 

o 


( I Oh) 


where H^,( D is the Planck flux (i e.. ?r times the usual Planck function). u t ,(p. p’i is the 
optical thickness defined by (61. and f (ui is the diffuse transmittance, wh h is proportional 
to the third exponential in.egral. 


flu) = 2 I- , <u). 


I n (.0 



(Ill 


In ( 10). p s denotes surface pressure, and l s denotes temperature ol the surface itself, which 
could differ from the air temperature immediately above the surface. T(p s ). 

The cooling rate ..t any level is related to the divergence of the net llu\ at that level 
integrated over wave number. 


01 ( pi 


( 12 ) 


Ot 





where C p is the specific heat of air at constant pressure and t is time. In the numerical 
examples considered here, the integration over wave number is carried out only for the por- 
tions of the water vapor spectrum listed in Table I . 

For computational convenience, we integrate ( lOai and ( I Ob) by parts and subtract one 
from the other to obtain 

Fj (P>- F„*(p) * B„(T<p t ))j f (u^fp. p,)) 

rT(p) dB,.(T) 

+ B„(T(( >) f (u„ (O.pi) + I r (u,, (p'.pi) — — — d I (p ) 

J T(0) T(p') 

where T(0i represents the temperature at the top of the atmosphere Ip %0). 

Making use of the far wing approximation, as embodied in ( 7 >. and defining the 
function 

G(w.T) = f f (wk„(p r , T r )j B„(T) d»> , (14) 

the cooling rate becomes 





y 


I he function G can he rapidly computed via the k -distribution method (Arking and 


Grossman. 1^72). We divid** the spectrum into n intervals ■ *»j - i»j _ j , i ■ I n. 

and let 11,1 I t he a suitable mean of 11, (1 > in the i 1 * 1 interval. | Mere there is an implicit 
assumption that for each interval ^ ( <»'<»',, JH, (T) H,, ( l>| and its derivative are 

smaller Ilian acceptable upper bounds | Letting x = log k. x,""" and Xj n,ax the minimum 
and maximum values of x in the i 1 ’ 1 interval, and hj (x) the distribution of x = log k (p r .T r ) 
normali/ed so (hat 


I. 


hj(x)dx* I. 


then the mean dljfusv transmittance in the interval, as a function of the scaled water vapor 
amount, becomes 


r, (wi 




Kj < wkt hj < x > dx 


( lb) 


fhe function (i can then be expressed in terms of r( w) as follows: 


n 

G(w.T) - Bj(T) f| (w) Ai^ 


(17) 


1 = 1 


Using the line parameters compiled by McOatchey el d. ( I ‘>72) and assuming a Voigt 
line profile, the k -distribution function, h(x), has been computed, and the results arc listed 
in Tables 3 and 4, respectively, for the band center am. the band wing regions at selected 
reference conditions. For repeated rapid computations one can store either the function 
fj ( w) or G ( w. T) and use it in a table lookup procedure. For the examples considered here, 
we employed the latter scheme. G is tabulated separately for the band center and band wing 
regions in Tables 5 and b. respectively. Fhe function behaves rather smoothly and is capable 
of providing high accuracy with minimal tabulation. 
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4. Comparison of methods 

a. Testing a p prt xi matio ns lin o by lino 

As a standard for testing various approximations, transmittances and cooling rales were 
calculated line by line and averaged over the appropriate spectral intervals. There are many 
ways of carrying out these calculations: line shapes and functional dependences have to be 
specified and appropriate numerical schemes must be* employed. Nevertheless, the line by 
line calculations employed here will be treated as “exact" for purposes of comparing the 
accuracies of various approximations and computational methods. The method for the line 
by line calculations will be fully described ! y Arking and Chesters ( l‘>7‘>) and briefly de- 
scribed as follows. 

The calculations are based upon the line parameters for the 11,0 molecule compiled by 
McClatchey et al. ( l‘>72). The line shape is assumed to be a Voigt profile cut off at 10 cm' 1 
from the line center, and the Lorentz half-width is assumed to be proportional to p /V' f . 

The atmosphere is assumed to consist of 30 layers overlying a blackbody surface, Each layer 
is treated as if it were isothermal at a temperature corresponding to the arithmetic mean of 
the temperatures at the upper and lower boundaries of the layer. The surface is assumed to 
be at 1000 mb and the top of the atmosphere at 0.2 mb. Using an accurate algorithm for the 
pressure integral of the Voigt 'unction developed by Milman ( I ‘>701. the integral ol the ab- 
sorption coefficient between each t air of levels is calculated monochromatically at intervals 
of .01 cm" 1 These integrals constitute the optical thickness between levels defined by tot. 
Monochromatic transmitance is then given by (5) and the cooling rates by ( 1 0 > and ( 12). 

The cooling rates were calculated line by line exactly as just described and with three 
different approximations: ( I > the far wing approximation, where Eq. (4) is used: ( 2) the 
isothermal approximation, where the absorption coefficient is computed at 250 k but pres- 
sure dependence is taken into account exactly; and (3) the diffusivity approximation . where 
the diffuse transmittance. 2 E^ (u). is replaced by exp (-1.66 u). 


Kig. 5a shows flu* exact cooling rates for the 200-580 cm -1 region. The lower wave 
numbers are near the center of the water vapor rotational band, an I the higher wave num- 
bers are in the wing of the band. It can be seen that the level of maximum cooling rate 
shifts upward as absorption increases. The differences between the exact cooling rate and 
that computed using the thiee approximations are shown in I'ig. 5b (approximation minus 
exact). It can bo seen that the far-wing approximation introduces only small errors with a 
maximum magnitude less than 2'? of the peak cooling rate. It is comparable to the error 
associated with the commonly used diffusivity approximation. In the 200-580 cm' 1 spectral 
region, the mean intensity transmittance, r. decreases with increasing water vapor amount 
w nearly according to the strong line approximation 

rtew ) * exp | -(cw) 1 '- | . 

and the mean diffuse transmittance, f. is approximately given by 

f(cw) - - / exp |-(cw//a)' J /i du| = 4 I ij |(ew)' /2 |, 

-X) 

where e is a function of temperature, pressure, and spectral interval. It is found that the 
difference between exp |-( 1 .60 cw )' /J | and 41 s l(cw) ,/2 | is less than 0.01 for all values 
of cw. Therefore, the use of the diffusivity approximation which amplifies the water vapor 
amount by a factor of ! .66 for the diffuse transmittance computations induces only small 
errors in the computed cooling rates. 

l or the isothermal approximation, the emission is correctly computed near the 400-mb 
level where the temperature is equal to 250 K for the tropical profile used. Since the tropo- 
spheric temperature decreases with increasing height, and since absorption decreases with de- 
creasing temperature, the atmosphere is computed too transparent below and too opaque above 
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(he 400-mb level. Accordingly, hoi li the computed incoming upward and downward fluxes 
at the 40()-mb level are larger than the exact ones, and more energy is absorbed near the 
400-mb level. The effect of over-estimatirg the energy absorbed near 400 mb reduces the 
cooling rate, as seen in Fig. 5b. The maximum error is I l i of the peak cooling rate, 
b. Present metho d compared with the Rodgcrs-Walshaw method 
l lte cooling rate calculations of Rodgers and Walshaw ( I %h) are based upon Cioody’s 
statistical band model for water vapor transmittance. The mhomogencity of the atmosphere 
is taken into account by using the Curtis-Godson approximation, the dilTusivitv approxima- 
tion is used for computing diffuse transmittances, and the spectrum is divided into a number 
of intervals, depending upon the desired accuracy. In this study we used 34 intervals, of 
40 cm*' each. (For application to atmospheric models, such as in the ( • I SS fK'M, the use of 
such fine intervals is considered too time consuming; instead, the intervals are 200 cm .1 
In the present method there is no need lor a large number of intervals. Instead, each 
portion of the spectrum is assigned to either of two region types, hand anii rs or hand wings 
and the spectral regions of each t\ pe are combined. The function G ( w. T) defined by 1 1 4) 
is pre-computed for each region type and the cooling rates are computed from (15). 

In order to compare the Rodgers- Walshaw method with the present method, the cooi.ng 
rate is computed exactly for the tropical and sub-arctic model atmospheres and is shown in 
Fig. 6a. It can be seen that the major contribution of the water vapor bands is in the mid- 
troposphere with maximum near the 300-mb level for the tropical model and 550 mb for 
the sub-a r ctic winter model. Near the surface and in the upper troposphere and above, the 
atmosphere is either too opaque or too transparent to have significant cooling. 

The errors (measured against the line by line calculations) in the cooling rate computed 
from the present method and the Rodgers-Walshaw method are shown in I igs. (>h and be for 
the tropical and sub-arctic winter cases, respectively The cooling rate error of the present 
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method is very small with a maximum magnitude of 4 V< ol the peak cooling rate, while the 
maximum error of the Kodgers-Walshaw method is 1 1'? of the peak cooling rate. 

The upward fluxes at the top of the atmosphere. I t (0), and the downward fluxes at 
the earth’s surface. I i (p N >. for the tropical and sub-arctic winter atmospheres are shown in 
Table 7. It can lx 1 seen that in both the present method and the Kodgers-W.dshaw method 
the errors are small. Because of the use of (4) in the present method, transmittances are 
computed too large in the upper troposphere and too small near the surface, and conse- 
quently the effective emitting levels for F t (0) and F i ( p s ) arc shifted downward. There- 
fore. both F t (0) and F J lp s ) are larger than the exact ones. For the Rodgers-Walshaw 
method, the computed transmittances are generally too small, and hence the value of Ft (0) 
is too small and that of I 1 (p t ) is too large compared to the exact ones, 
c I he effect of the e-typ e continuous absor pt inn 

The results shown in this study were obtained by including only the absorption by 
water vapor lines; e-type continuous absorption was not considered. We have investigated 
the effect of e-type absorption in the 400-5K0 cm" 1 spectral region lor the tropical model 
atmosphere. Using the line by line calculation, it has been found that the effect is to reduce 
the transmittance and. hence, the cooling rate by 0.1 5° Cl day in the ‘>00-1000 mb region. 
Above the 900-nib level, the effect is very small. For drier atmospheres, the effect of the 
e-type absorption in the water vapor bands can be neglected. In another study (Chow and 
Arking. I l >78). e-type continuum absorption and the 10 pm window region were included 
in the computation of cooling rates. 
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5. Summary mul further remarks 

An accurate hut fast method has been developed for infrared llux computations in the 
water vjpor hands, file absorption coefficient at any part of the atmosphere is computed 
from the absorption coefficient at the reference conditions using the far-winy* approximation 
lor which the dependence of the absorption coefficient on pressure and temperature is as- 
sumed to correspond to that of the wings of absorption lines. The present method differs 
from the commonly used one-parameter scaling by suitably choosing the reference conditions 
corresponding to the regions where cooling is most significant. With the use ol the far- wing 
approximation, the wavenumber dependence of the absorption coefficient can be separated 
from pressure and temperature, and the computation is greatly simplified. Further reduction 
of the computation is achieved by dividing the entire water vapor spectrum into only two 
regions (band center and wing regions) and using the pre-corn pu ted tables. 

Using the exact line-by-line calculation as a standard, the accuracy of the present 
method is compared with that ol Rodgers and Walshaw. It has been found that the present 
method introduces a maximum error of only 4 'a ol the peak cooling rate, while the Rodgers- 
Walshaw method introduces an W'Z error. The upward (luxes at the top of the atmosphere 
and the downward fluxes at the earth’s surface were also compared with the exact line-by- 
line calculations; the errors in these fluxes are small with a magnitude less than 2 '< for both 
the present method and the Rodgers Walshaw method. 

The accuracy of the far-wing approximation was also compared with the commonly 
used diffusivity and isothermal approximations. I or the diffusivity approximation, the 
diffuse transmittance was computed as the intensity transmittance but the optical thickness 
was amplified by a factor of 1.66. For the isothermal approximation, the transmittance 
was computed at a temperature of 250 K it has been found that both the far-wing and the 
diffusivity approximations introduce only small errors in the oinputed cooling rates with 


15 



a maximum error !$ 2 !% ol the peak cooling rale, while the isothermal approximation intro- 
duces large errors ( ~ 1 1 ). The cooling rate for the isothermal approximation is under- 
estimated in the neighborhood of the 250 K level and is overestimated near the surface and 
in the upper troposphere. 
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Tabic I. The reference temperature T f and pressure p f are specified for the center and the 
wing regions of the water vapor absorption bands. The effect of temperature on the 
absorption coefficient is indicated by K (T. T f ). which is the mean value of the function 


R,,(T.T r ) defined by l.q. (3). 



Wave number (cm* 1 ) T f < K ) 

p r (mb) 

R(T r - 40 K. T r ) 

R«T r ♦ 40 K. T r ) 

Center 

0-340 'l 
1380-1000 J 

► 225 

275 

0.00 

1.16 


140-580 ^ 





Wing 

1220- 1380 

► 25b 

550 

0.58 

1.78 


1000-2020 j 







Table 2. Mean 

intensity transmittance over 20 cm 

spectral intervals. 


l.aycr 

(mb) 

320-340 cm” 

"Exact" 

transmittance 

1 

Model 

transmittance 

Layer 

(mb) 

460-480 cm -1 

"Exact” Model 

transmittance transmittance 

10-100 

0.010 

0.015 

10100 

0.987 

0.996 

100-250 

0.814 

0.815 

100-300 

0.95 1 

0.957 

250-400 

0.350 

0.357 

300-600 

0.530 

0.531 

400-500 

0.1 16 

0.125 

600-850 

0.162 

0.161 

500-600 

0.025 

0.036 

850-1000 

0.100 

0.099 


10 



fj 

c 


Table 3. The k-distrihuticn function for the band center region with p r = 275 mb and T f = 225 K. 


Wave n jmber 


20- 60 

0.0 

C.O 

0. 0 

0. 0 

0. 0 

0. 0 

60- ICC 

0. 0 

o.c 

0. c 

C.O 

0.0 

0. 0 

ICC- 190 

0.0 

c.c 

o.c 

0.0 

0. 0 

0. 0 

19C- 18C 

0.0 

o.c 

c. 0 

0.0 

0. 0 

0. 0 

7 SC- 220 

C.O 

O.c 

o.c 

0.0 

0.0 

0.0 

220- 260 

c.c 

0.0 

0.0 

0. 0 

0. 0 

0.0 

260- 3C0 

C.O 

o.c 

0.0 

0.0 

0.0 

0.0 

3C0- 390 

o.c 

0. 0 

0.0 

0.0 

C.O 

0.0 

•l 350- 1 92C 

c.c 

C.O 

0.09 1 

C. 288 

0.933 

0.577 

•1920-1960 

O.C 

0.0 

o.c 

o.c 

0. 203 

0. 309 

1196 0-15CC 

C.O 

0.0 

0. 0 

0.0 

0.0 

0. 010 

15CO-159C 

o.c 

c.c 

o.c 

C.O 

0. 0 

C.O 

1590-1560 

c.c 

o.c 

o.c 

c.c 

0.0 

c. 0 

158C- 1670 

0.0 

0.076 

0. C?9 

0. 155 

0.210 

0.69C 

7620-1660 

o.c 

c.c 

0. 0 

O.C 

0.0 

0.0 

1660-1700 

0. 0 

C.O 

0. 0 

0.0 

0.0 

0. 0 

1 700- 1790 

c. 0 

0. 0 

0.0 

0.0 

0.0 

0.0 

179C-178C 

o.c 

0. 0 

0.0 

0. 0 

0. 0 

0.0 

1780- 182C 

o.c 

c.c 

0.079 

0. 362 

0. 161 

0.975 

182C-1P6C 

0. 1 19 

0. 108 

0.098 

0.319 

0. 297 

0.529 

1860-1900 

0.086 

0.395 

0. 386 

0.679 

0.669 

0. 997 


lo <»10 k 

0.50 0.75 1.00 1.25 1.50 I.75 

0.101 0.«l76 0. 263 0. 505 0.597 0. 301 
0.0 0.0 0.0 0.0 0. icg 0.2511 

0.0 0.0 0.0 0.0 0.253 0. 102 

0.0 0.0 C.O O.c 0.331 0.919 

0. 0 0. 0 0. 0 0. 097 0.256 0. HO 

C.0C8 0.379 0. 379 0.301 0.223 0.268 
0.0 0.686 0.707 0.593 0.969 0.239 

0.300 0.798 0.661 C.930 O.933 O.335 
0.585 0.999 0.337 0.339 0.259 0.170 
0.927 0.662 0.653 0.38; 0.279 0.226 
0.339 0.366 0.721 0.607 0.966 0. 395 
0.0 0.0 0.169 0.932 0.782 0.607 

C.O C.O 0.272 0.722 0.602 0.991 
0.999 0.603 0.928 0.331 0.276 0.209 
C.O 0.0 0.570 0.757 0.959 0.528 

0,0 0-0 0.299 0.680 O.638 0.638 

0.178 0.320 0.952 O.538 0.966 0.506 
0.162 0.695 0.508 0.759 0.932 0. 39C 
0.975 O.990 0.501 0.385 0.275 0.217 
0.527 0.557 0.370 0.269 0.196 0.199 
0.279 0.179 0.220 0. 160 0. 193 0. 103 


2.00 
0.230 
0.513 
0. 969 
0.901 
0. 199 
0.935 
0.212 
0 . 221 
0. 127 
0.227 
0. 329 

0. 355 
0. 920 

0. ne 
0 . 386 
0. 998 
0. 356 
0.290 
0. 167 
0. 197 
0. 070 


2.25 
0. 321 
0.682 
0.815 
0.537 
0. 357 
0.501 
0. 355 
0. 159 
0.088 
0. 185 
0.220 
0. 385 
0. 352 
0. 1 13 
0. 322 
0.332 
0.287 
0.226 
0. 115 
0.089 
0.051 


2.50 
0.333 
0.602 
0.511 
0.512 
0.261 
C.397 
0.227 
0. 1J6 
0.088 
0. 103 

0. 1t>2 

C. 323 

0.269 
o.oec 
0.216 
C. 257 
0.236 
0. 170 
0.079 
0.069 
0.038 


Table 3. (Continued) 



Wave number 

(cm' • ) 

2.75 3. CO 3 

1.25 3 

1.50 3 

.75 « 

.00 “ 

: °9l0 h 

.25 “.50 “ 

.75 5 

.00 5 

.25 5 

.50 5 

.75 6.00 6 

.25 

20- 6 C 

C. 213 

0. 191 

C. 136 

C. 101 

0.083 

C.C73 

O.C53 

C.C3C 

0.024 

0.01 1 

C. 006 

O.C 

C.C O.C 

c.c 

to- IOC 

C.“63 

0. 389 

0. 303 

0. 172 

0. 128 

0. 1 19 

0. 088 

0.069 

0.0*3 

0.039 

0.01“ 

0.CC8 

C.CC5 C.O 

C.C 

ico- me 

C. 37? 

0. 32fc 

C. 322 

C. 220 

0. 167 

o. i 2 e 

0. 1C7 

0. 072 

C.C“9 

C.C35 

C.C29 

C.C16 

C.C12 C.C 

O.C 

me- isc 

C. 1IH8 

C.351 

0.25 3 

0. 177 

0. 176 

0. 122 

c.ceo 

O.G59 

0. 0“ 3 

C.C“C 

0.02“ 

0.C18 

C.CC7 C.CC6 

C.C 

18C- 22 0 

C.2C7 

0. 156 

0. 097 

0.070 

0.052 

O.Oi* 

C. 026 

0.C1S 

0. 026 

0.C17 

0.01“ 

C. C 1 3 

C.CG5 O.CC* 

0* 0 

- 

C.?t7 

0. 230 

0. 172 

C. 125 

C.C90 

0.061 

0. 0“3 

C.C3* 

0.035 

0.023 

C.015 

e.cco 

C.Cf'? C.CC3 

C.C 

• - • 

C. 15“ 

C. 1C5 

0.069 

C. 063 

0. 0“2 

0.0*1 

C.C3C 

0.02“ 

C.C15 

0.012 

C.CC* 

C.CC? 

c.c c.c 

c.c 

- 3“C 

C. It 1 

C. 09“ 

C. 063 

C. 063 

0.0*1 

C. 027 

0. 026 

C. 019 

0.020 

C.CC7 

0.C03 

C.CC3 

c.c c.c 

c.c 


O.C65 

0.C“1 

C.C35 

C. 029 

0.023 

C.C1“ 

0. CCS 

. 

0.0 

c.c 

C.C 


C.C c.' 

c.c 

i* 2 c-i“ 6 o 

C. 075 

C.C62 

0.066 

0.C60 

C.03C 

C.C26 

C.C1C 

0. 008 

o.oc* 

O.C 

U.C 

O.C 

C.C c.c 

C.C 

• • • - • 

C. 121 

C.09C 

G. 065 

0.C51 

0. 036 

C. 032 

C. 027 

0.016 

C.CC2 

c.c 

O.C 

C.C 

c.c c.c 

rs 

15CC-15“C 

C. 3C 3 

0. 19“ 

0.139 

0. 10“ 

0. 076 

0.055 

0.03“ 

0.02* 

0.017 

0.CC6 

O.C 

c.c 

c.c c.c 

C.C 

15“C-158C 

0.22“ 

C. 177 

C. 121 

C. ICO 

0.07“ 

C. 063 

C.053 

0.C28 

0.020 

C.C12 

n.O 

C.C 

C.C O.C 

c.c 

1580-1620 

C.C65 

e*“ 

0.031 

C. C 1 5 

o.cn 

C.C06 

0. 00c 

0.00“ 

c.cc* 

- 

C.C 

O.C 

c.c c.c 

c.c 

162C-166C 

C. 177 

0. 155 

0. 125 

C. 066 

0. 058 

C. 052 

C. 0“6 

0.C29 

C.C15 

0.C12 

O.C 

c.c 

c.c c.c 

c.c 

1 € € 0 - 1 7C0 

C. 1 E 2 

c. i?e 

0. 128 

0. 107 

0.C71 

C.0“8 

C.C*C 

0.02“ 

0.C15 

0.C12 

C.CC2 

c.c 

c.c c.c 

c.c 

-17H 

C. 1 3? 

0. 1“3 

0. 101 

C. 087 

0. 075 

0. 039 

0.038 

0.C26 

C.C1? 

C. 002 

C.C 

c.o 

c.c c.c 

c.c 

17*0-1730 

C. 123 

C.C87 

0. 072 

0. C“ 3 

0.038 

0.021 

C.C20 

0.010 

C.CC« 

O.C 

O.C 

O.C 

c.c c.c 

c.c 

• • . ■ 

C.C66 

C. C“9 

0.038 

0. 029 

0.C2C 

0.C12 

O.CC« 

c.co* 

C.CC 2 

C.C 

c.c 

c.c 

c.c c.c 

C.C 

-18( 

C. C“7 

C.C3’ 

C.C“C 

C. 017 

0.018 

0.01“ 

C.C06 

0.00“ 

O.C 

C.C 

2.0 

c.c 

c.c c.c 

c.c 

' 

C.C35 

0.027 

0. C 1 5 

0.012 

C. 008 

0.0C6 

0.0C2 

c.o 

0.0 

C.O 

O.C 

. 

O.C O.C 

O.C 



Tabic 4. The k -distribution function for the hand *mg region with p r = 550 mb and T ( 


256 k 


Wave number 

(em'M 

-2.50 

-2.25 

-2.00 

-1.75 

-1.50 

-1.25 

o 

o 

a 

IO9, 0 k 

-0.75 ■ 

-0.50 ■ 

-0.25 

0.0 

0.25 

0.50 

0.75 

1.00 

3*0- 390 

O.C 

C.O 

O.C 

C.O 

0.0 

O.C 

O.C 

0.0 

0.0 

O.C 

O.C 

0.26C 

C. 299 

0.2*2 

C. 569 

380- «20 

0.0 

c.c 

0.0 

O.C 

O.C 

O.C 

O.C 

C.C70 

0. to* 

C.261 

0.335 

0. *80 

0.623 

0.399 

0.2*3 

*2C- *60 

0.0 

0.0 

0.0 

c.o 

0.0 

C.O 

C.O 

0.0 

O.C 

C.08C 

0.525 

0.672 

C.51 1 

0. *98 

0. *93 

*ec- 500 

O.C 

c.o 

c.c 

0.0 

0.0 

0.0 

c.o 

0.231 

0.297 

0.815 

0.531 

0.583 

C. *77 

C.3C0 

3. 19* 

5C0- 5*0 

0.0 

0.0 

0.0 

0.0 

C.C 

0.0 

C. 123 

c. 380 

0.369 

0.5*7 

0. *72 

0.5** 

C. *07 

0 . 32 * 

0.205 

5*C- 580 

0.0 

0.0 

c.o 

0.0 

0.0 

C. 1 15 

C. 199 

0. *08 

0.856 

0.602 

0. *3* 

0.*C3 

0.2*7 

0.2C5 

0. 15* 

1220- 1260 

0. 196 

0.216 

C 12C 

0. *70 

0.605 

0.609 

0. #89 

0. 369 

0.280 

0.227 

0. 169 

0. 1 12 

0.055 

0. 030 

0.018 

1260-13CC 

0.0 

c.o 

C.O 

0.0 

0.212 

0. *76 

0.51* 

0.527 

o.* e * 

0. * 1 1 

0 . *32 

0.276 

0.203 

0. 162 

0.111 

1300-13*0 

0.0 

c.o 

C.O 

0.09* 

0. 152 

C. 121 

C.C79 

0.200 

0. 368 

0.588 

0.517 

0.326 

0. 36e 

0.336 

0.212 

1 3*C- 1 380 

0.0 

c.o 

c.o 

0.0 

0.0 

C.C 

0.0 

0.0 

0.070 

0.361 

0. *81 

C. *73 

0.705 

0. **3 

0.332 

1 9CC- 1 9* 0 

O.C 

0.0 

O.C 

c.o 

0.0 

0.0 

0.0 

C.5S» 

C. 2*2 

0. *2C 

C.739 

c. * 17 

0. *50 

0.327 

0.220 

19*C-1?5o 

0.0 

0.0 

0.0 

0. 103 

0.082 

0. 116 

0.336 

9.500 

0.562 

0. *70 

C. 362 

C.289 

0.277 

0. 197 

0.205 

19 S 0 - 2 C 2 C 

c.o 

0.0 

O.C 

0.0 

0. 183 

C.*38 

0.579 

C .533 

0.507 

0.322 

0.330 

0.288 

0.20* 

0. 15* 

0. 107 


Table 4. (Continued) 


-4 'J 


Wave number 
(cm* 1 ) 

3 « 0 - 380 
380- 920 
920- 60 

960- 500 
50C- 5^0 
59C- 580 
1220-1260 
1260-1300 
1 300- 1 3*‘0 
13HC-1330 
1900-1990 
1 9 J 4 0 - 4 950 
1980-2020 


1.25 

1.50 

1.75 

2 . 0 C 

2.25 

2.50 

2.75 

*°9l0 k 

3.00 3.25 

3-50 

3.75 

9 . 0 c 

9.25 

5.50 

9.75 

C.996 

0.263 

0. 302 

C. 91|A 

0. 359 

0.213 

0. 137 

0. 119 

0.067 

0.073 

0.061 

C.C96 

0.023 

0.012 

C.CC 

C. 197 

0.2C9 

C.218 

0 . 135 

0.216 

0. 199 

c .111 

0.083 

0.C60 

0.099 

0.039 

C. 022 

C. 007 

O.C 

C.C 

0. 327 

0.205 

c. 19 c 

0. 129 

0. 1 12 

0.C80 

0.063 

O.C63 

0.028 

0.011 

0.005 

0.CC3 

O.C 

O.C 

C.C 

0. 192 

0. 1 15 

0.C85 

0.C51 

C.C59 

0.051 

O.C33 

0.018 

0.017 

0.006 

0.0 

C.O 

C.C 

C.O 

O.C 

C. 188 

c. 1 30 

0.09C 

0.C65 

0.09*i 

C.C 36 

0.026 

0.029 

0.013 

0.008 

O.C 

O.C 

C.C 

O.C 

C.C 

0. ICO 

c.ceo 

C.C60 

C. 05*1 

C. 035 

0.030 

0.C13 

C. 003 

0.002 

O.C 

O.C 

C.O 

C.C 

O.C 

C.O 

C.C15 

C.C13 

c.cce 

C.C 

C.C 

O.C 

C.C 

O.C 

0.0 

C.C 

0.0 

C.O 

C.C 

O.C 

C.C 

c.o89 

C.C59 

o.cuc 

0.C26 

0.012 

C.0C6 

C.C 

C • 0 

0.0 

0.0 

O.C 

C.O 

C.C 

C.C 

C.C 

0. 137 

0. 165 

0. 115 

C . 079 

C.C62 

0.039 

0.021 

0.016 

C.0C3 

O.C 

C.C 

C.O 

C.O 

0.0 

C.C 

C . 295 

C.212 

A . 159 

0. 1 35 

0. 101 

C.079 

C. 06C 

0.096 

C.033 

0.016 

0.009 

C.C 

O.C 

C.O 

0.0 

c. 132 

C. 138 

0.C96 

0.072 

C.C**8 

C.C37 

0.03C 

C. 029 

0.C1 1 

0.009 

0.CC9 

*0.005 

C.O 

0.0 

C.C 

C. 1 tic 

U. 1C? 

0.C69 

C.C55 

C.C91 

0.C37 

C. 029 

0.016 

0.006 

0.006 

0.0 

C.C 

C.C 

0.0 

C.C 

c.eec 

0.C97 

0.058 

0.052 

C.C28 

0.011 

0.013 

0.005 

0.009 

0.002 

0.0 

O.C 

O.C 

O.C 

C.C 



Table 5. Values of <i(w, T,p r , 1 f ) lor the band center region, in which p f and T f are set, 
respectively, at 275 mb and 225 K. The units are C;<W/m : ), T(K). w(g/cm 2 ). 


|o g| ( ,W 


200 


210 


Temperature (K) 

220 230 240 250 


-5.0 

24.77 

27.89 

-4 8 

24.51 


-4.6 

24.20 

27.25 

-4.4 

23 80 

26.82 

-4.2 

23.32 

26.28 

-4.0 

22.73 

25.63 

-3.8 

22.02 

24.84 

-3.6 

21.17 

23.89 

-3.4 

20.15 

22.77 

-3.2 

18.96 

21.45 

-3.0 

1 7.60 

19.94 

-2.8 

16.08 

18.24 

-2.6 

14.42 

16.40 

-2.4 

12.68 

14.45 

_ s ■> 

10.91 

12.47 

-2.0 

9.17 

10.50 

-1.8 

7.49 

8.61 

-1.6 

5.92 

6.82 

-1.4 

4.46 

5.16 

-1.2 

3.15 

3.67 

-1.0 

2.04 

2.39 

-0.8 

1.17 

1.39 

-0.6 

0.57 

0.70 

-0.4 

0.23 

0.30 

-0.2 

0.08 

0.12 

-0.0 

0.03 

0.05 

0.2 

0.01 

0.03 

0.4 

0.01 

0.01 

0.6 

0.00 

0.01 


31.36 

35.29 

39.79 

31.05 

34.95 

39.41 

30.66 

34.52 

38.94 

30.18 

33.99 

38.35 

29.59 

33.34 

37.64 

28.88 

32.55 

36.77 

28.00 

31.59 

35.70 

26.95 

30.43 

34.42 

25.71 

29.05 

32.90 

24.25 

27.44 

31.12 

22.57 

25.58 

29.06 

20.69 

23.49 

26.75 

18.64 

21.21 

24.21 

16.47 

18.79 

21.51 

14.24 

16.31 

18.73 

12.04 

13.83 

15.95 

9.91 

1 1.42 

13.22 

7.88 

9.12 

10.60 

5.99 

6.97 

8.14 

4.28 

5.01 

5.90 

2.82 

3.34 

3.97 

1.67 

2.01 

2.44 

0.86 

1.08 

1.35 

0.39 

0.52 

0.68 

0.17 

0.24 

0.34 

0.08 

0.12 

0.18 

0.04 

0.06 

0.10 

0.02 

0.03 

0.05 

0.01 

0.02 

0.02 


44.97 

44.56 
44.04 

43.40 
42.61 
41.65 
40.47 
39.06 

37.37 

35.40 
33.12 
30.55 
27.73 
24.71 
21.59 
18.45 
15.36 

12.38 

9.56 

6.98 
4.76 

2.98 
1.70 
0.91 
0.48 
0.26 
0.14 
0.08 
0.04 
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Tabic 5. (Continued) 


OF /), , , ' /•; j 

° li WAUT Y 


Temperature (K) 


-5.0 

50.99 

58.01 

-4.8 

50.54 

57.51 

-4.6 

49.97 

56.88 

-4.4 

49.26 

56.10 

-4.2 

48.39 

55.14 

-4.0 

47.33 

53.96 

-3.8 

46.03 

52.52 

-3.6 

44.47 

50.79 

-3.4 

42.60 

48.72 

-3.2 

40.41 

46.29 

-3.0 

37.88 

43.48 

- 2.8 

35.02 

40.29 

- 2.6 

31.88 

36.77 

-2.4 

28.50 

32.98 

- 2.2 

24.99 

29.01 

- 2.0 

21.43 

24. 9K 

- 1.8 

17.91 

20.97 

- 1.6 

14.50 

17.05 

-14 

11.26 

13.32 

- 1.2 

8.30 

9.89 


6 ( 1 . 1 9 
65.63 
64.94 
64.07 

63.01 
61 70 

60. 1 1 
58. IS 
55.88 

53.18 
50.04 
46.48 
42.53 
38.27 

33.78 

29.20 
24.60 

20.09 

15.78 
1 1.81 

8.34 
5.52 
3.42 
2.03 
1.19 
0.70 
0.40 
0.22 
0.1 1 


75.70 

75.09 

74.32 

73.36 

72.18 

70.72 
68.95 

66.81 
64.25 
61.23 

57.73 

53.74 
49.30 
44.49 
39.41 
34.17 
28.90 

23.71 

18.72 

14.10 
10.06 

6.76 

4.29 

2.61 


86.73 

86.05 

85.20 
84.14 

82.82 

81.20 
79.23 
76.84 
73.99 
70.61 
66.69 
62.22 
57.22 

51.78 
46.00 

40.02 

33.97 

27.97 

22.19 

16.83 
12.13 

8.26 

5.34 

3.32 

2.03 

1.22 

0.71 

0.39 

0.19 
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Table 6. 

Values of (i ( w 

. T p r . T 

r ) for the band wing region, in which p ( and T ( 

are set. 

respectively, at 550 

mb and 

256 K. The 

units are (i (W 

in* ). T(K), w (g/cm*). 




Temperature ( K 1 



IoKiqW 








200 

210 

220 

230 

240 

250 

-4.4 

33.13 

39.48 

46.51 

54.27 

62.80 

72.16 

-4.: 

32.99 

39.31 

46.32 

54.06 

62.56 

71.90 

-4.0 

32.80 

39.10 

46.07 

53.77 

62.24 

71.54 

-3.8 

32.55 

38.81 

45.75 

53.40 

61.82 

71.07 

-3.6 

32.24 

38.45 

45.33 

52.93 

61.29 

70.47 

-3.4 

31.85 

37.99 

44.80 

52.33 

60.61 

69.71 

-3.2 

31.36 

37.43 

44.15 

51.58 

59.77 

68.76 

-3.0 

30.77 

36.74 

43.36 

50.67 

58.74 

67.60 

-2.8 

30.06 

35.90 

42.39 

49.57 

57.49 

66.20 

-2.6 

29.21 

34.9 1 

41.25 

48.27 

56.01 

64.52 

-2.4 

28.22 

33.76 

39.9 1 

46.73 

54.26 

62.56 

_ T 2 

27.09 

32.43 

38.37 

44.97 

52.25 

60.29 

-2.0 

25.80 

30.92 

36.62 

42.95 

49.95 

57.68 

-1.8 

24.35 

29.21 

34.64 

40.67 

47.35 

54.73 

-1.6 

22.72 

27.29 

32 40 

38.09 

44.40 

51.38 

-1.4 

20.90 

25.14 

29.89 

35.18 

41.07 

47.59 

-1.2 

18.87 

22.74 

27.08 

3 1 .94 

37.34 

43.34 

-1.0 

16.67 

20.13 

24.02 

28.38 

33.26 

38.68 

-0.8 

14.33 

17.35 

20.76 

24.60 

28.90 

33.69 

-0.6 

1 1 .94 

14.51 

17.41 

20.69 

24.39 

28.52 

-0.4 

9.59 

1 1 .69 

14.09 

16.82 

19.90 

23.36 

-0.2 

7.38 

9.04 

10.95 

13.14 

15.62 

18.44 

0.0 

5.41 

6.67 

8.14 

9.83 

1 1.76 

13.97 

0.2 

3.76 

4.68 

5.76 

7.01 

8.47 

10 14 

0.4 

2.47 

3.10 

3.86 

4.76 

5.81 

7.04 

0.6 

1 51 

1.93 

2.44 

3.06 

3.80 

4.67 

0.8 

0.87 

1.13 

1.47 

1.88 

2.38 

2.97 

1.0 

0.47 

0.63 

0.85 

i . ! 1 

1.44 

1.85 

1.2 

0.25 

0.35 

0.48 

0.65 

0.87 

1.14 

1.4 

0.13 

0.20 

0.28 

0.39 

0.53 

0.71 


Table 6 . ( Continued ) 


i 'AGE IS 

• 1 m QUALITY 


log I0 w 

260 

270 

280 

290 

300 

310 

- 4.4 

82.41 

93.60 

105.80 

1 19.08 

133.50 

149.14 

- 4.2 

82.11 

93.27 

105.44 

1 18.69 

133.08 

148.67 

- 4.0 

81.71 

92.83 


118.16 

1 32 50 

148.04 

- 3.8 

81.19 

92.26 

104.32 

1 1 7.46 

131.74 

147.22 

- 3.6 

80.52 

91.51 


1 16.56 

130.75 

146.14 

- 3.4 

79.67 

90.57 

102.47 

1 15.42 

1 29.50 

144.78 

- 3.2 

78.62 

89.40 

101.17 

1 13.99 

127 94 

143.07 

- 3.0 

77.32 

87.96 

99.58 

1 12.24 

126.01 

140.96 

1 

1 J 

be 

75.75 

86.2 1 

97.64 

1 10.10 

123.66 

138.39 

- 2.6 

73.88 

84.12 

95.32 

107.55 


135.31 

- 2.4 

71.67 

8 1 .66 

92.60 

104.53 

1 17.53 

131.66 


69.12 

78.81 

89.43 

101.02 

1 13.66 

127.41 

- 2.0 

66.19 

75.54 

85.78 

96.98 

109.19 

122.49 

- 1.8 

62.86 

71.81 

8 1 .62 

92.35 

104.07 

1 16.84 

- 1.6 

59.08 

67.56 

76 87 

87.07 

98.22 

1 10.38 

- 1.4 

54.79 

62.74 

71.48 

81.07 

91.56 

103.01 

- 1.2 

49.99 

57.33 

65.42 

74.30 

84.04 

94.68 

- 1.0 

44.70 

51.37 

58.72 

66.82 

75.71 

85.45 

- 0.8 

39.03 

44.96 

51.52 

58.76 

66 73 

75.46 

- 0.6 

33. 1 5 

38.30 

44.01 

50.34 

57.32 

64.99 

- 0.4 

27.26 

3 1 .6 1 

36.46 

4 1 .86 

47.82 

54.40 

- 0.2 

21.62 

25.19 

29.20 

33.66 

38.63 

44 12 

00 

16.49 

19.33 

22.54 

26.14 

30.15 

34.61 

0.2 

12.07 

14.20 

16.75 

19.57 

i ■> -j ■) 

26.24 

0.4 

8.46 

10.10 

1 1.98 

14.1 1 

16.52 

19.22 

0.6 

5.69 

6.88 

8.25 

9.82 

1 1.61 

13.63 

0.8 

3.69 

4.53 

5.51 

6.64 

7.94 

9.41 

1.0 

2.34 

2.92 

3.60 

4.39 

5.31 

6.35 

1.2 

1.47 

1.86 

2.33 

2.87 

3.50 

4.21 

1.4 

0.92 

1.18 

1.49 

1.85 

2.27 

2.75 









Table 7. Upward fluxes at the top of the atmosphere, F t (01. and downward fluxes at the 
earth’s surface. F J ( p s ), in W/m*. Percentage errors in fluxes ate given in the parenthesis. 



Tropical 

Sub-arctic winter 


F t (0) 

F4(p s ) 

F t (0) 

n<p s > 

Exact 

1 33 4U 

21 1 .88 

106.85 

99.45 

Present method 

135.42 

213.33 

108.09 

100.39 


(+1.590 

( +0. 1 % ) 

( + 1.29?) 

1 + 1.090 

Rodgers-Walshaw method 

132.44 

2I5.‘)7 

106.61 

i00.38 


t-0.7'; > 

(+1.990 

(-0.290 

(+1 09?) 




10.0 


5.0 



Wave Number (cm 1 ) 


Figure lb. Effect of temperature on the absorption coefficient. R, (T, T, ), in the band wing region with T r = 256 k. The 
dashed curve is for T = T r + 40 K. and the solid curve is for T = T r - 40 K. 
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Figure 2a. Temperature profiles used in this study. The tropical and subarctic- 
winter profiles are from McClatchcy et al. (1972) 
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Pressure (mb) 
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Figure 2b. Humidity profiles associated with the temperature profiles in Figure 2a. 
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Figure 3a. Intensity transmittances between the 100- and 250-mb leveis for the spectral region 320-340 cnr> A :> p.cal 
miu -latitude atmosphere was used lor the computations. The solid curve was computed using the exact line-by- 
line method, and the dashed curve was computed using the far-wing approximation. 
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Figure 3b. 


Intensit\ transmittance* between the 250- and 400-mb levels for the spectral region 320-340 cnr» . A typical 
mid-latitud: atmosphere was used for the computations. The solid curve was computed using the exact line-hy 
line method, and the dashed curve was computed using the far-wing approximation 




Transmittance 



Wave Number (cm 1 } 


Int.nMty t admittances between the 300- and oOO-mb levels for the spectral regron 460-480 ciir » \ typical 

mid-latitude atmosphere was used for the computations. Die solid curve was computed using the exact hne-bv- 
•ne method, and the dashed curve was computed using the far-wing approximation * 
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Figure 5a. Cooling rates computed from the exact line-by-line method in the spectral 
region 200-580 cm"' . The clear tropical atmosphere shown in I igure 2 
Wu s used for the computations. 
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Figure 5b. Differences between the exact cooling rate and the cooling rates com- 
puted by using, respectively, the far-wing, the diffusivity, and the iso- 
thermal approximations for the spectral region 200-580 cm* 1 . 
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Figure 6a. Fxact cooling rates in the 0-580 cm* 1 and 1 220-2020 cm** spectral 
regions tor the clear tropical and subarctic-winter atmospheres shown 
in Figure 2. 
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Figure 6b. 


Differences between the exact cooling rate and the cooling rates com- 
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Figure be. Differences between the exact cooling rate and the cooling rates com- 
puted by using, respectively the present method and the Rodgets-Wal- 
shaw method in the 0-580 enr' and 1 220-2020 cm"' regions for the 
subarctic-winter atmosphere shown in Figure ? 
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